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Abstract 

We consider a model for a Z'-boson coupled only to B-L and hypercharge. Besides the usual right- 
handed neutrinos, we add a pair of fermions with a fractional lepton charge, which we therefore 
call leptinos. One of the leptinos is taken to be odd under an additional Z2 charge, the other even. 
This allows for a natural (inverse) seesaw mechanism for neutrino masses. The odd leptino is a 
candidate for dark matter, but has to be resonantly annihilated by the Z'-boson or the Higgs-boson 
responsible for giving mass to the former. Considering collider and cosmological bounds on the 
model, we find that the Z'-boson and/or the extra Higgs-boson can be seen at the LHC. With 
more pairs of leptinos leptogenesis is possible. 
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I. INTRODUCTION 



The standard model (SM) gives an excellent description of the known laws of particle 
physics. However there are a few facts that it cannot explain. Of importance for this 
paper are the neutrino mass spectrum, in particular the smallness of neutrino masses, the 
baryon-antibaryon asymmetry in the universe and the presence of dark matter. 

Though it appears to be very difficult to explain the mass-spectrum of the fermions, there 
is a mechanism that can in principle make it plausible why the neutrino masses are much 
smaller than the other masses. The mechanism uses the presence of right-handed neutrinos 
with a large Majorana mass. The mass-matrix also contains a Dirac-mass. If the Majorana 
mass is much larger than the Dirac mass, one finds after diagonalization one (very) light 
and one heavy neutrino. The method is called the seesaw mechanism. Several realizations 



exist, the simplest being the type-I seesaw 



5|. In this paper we consider a somewhat 



more mvo 
literature 



ved form, which is conventionally called the inverse seesaw mechanism in the 



Id 



The presence of heavy neutrinos does not only affect the mass-spectrum, it can also 



provide a mechanism to exp 
"leptogenesis" js] (also see 



ain the baryon asymmetry of the universe. The method is called 



12| for recent reviews). An asymmetry in the lepton sector is 



produced via CP-violating heavy lepton decays. Subsequently the asymmetry in the lepton 
sector is transfered to the baryons by means of electroweak (EW) sphaleron processes 
151] • Most early papers used neutrino masses at the grand unified scale, however a successful 
TeV scale leptogenesis is possible as well if two such neutrinos are degenerate in mass, 
thereby enhancing the CP-asymmetry parameter. This goes under the name of "resonant" 



leptogenesis 



16l-l2l|. In this case, fiavour effects have to be taken into account |22|-|26|. 



Cosmological observations imply the existence of non-baryonic dark matter that drives 
structure formation on large scales and dominates galactic and extra-galactic dynamics. 
Within the SM there is no candidate for this matter. One therefore has to enlarge the 
SM. The easiest way to explain the dark matter is to postulate the existence of thermally 
produced Weakly Interacting Massive Particles (WIMPs), with masses roughly at the weak 
scale. An additional unbroken symmetry is postulated, that prevents the WIMPs from 
decaying. 



The latest experimental searches (see, e.g., 27(|) confirm the SM again at higher energy 



scales than before. The data leave little space for modifications. In particular complicated 
extensions of the SM lead to phenomenological problems, for instance with flavour changing 
neutral currents, and tend to need many fine tunings of parameters. Minimal extensions 
are therefore preferable. The simplest form to enlarge the gauge group of the SM is to add 
a single U{1) factor, which has to be a linear combination of hypercharge and B-L, baryon 
minus lepton number, if one does not want to enlarge the fermion spectrum. These are 
the so-called "minimal Z' models" (see, for example 28|-|33|). The spontaneous symmetry 
breaking of the extra U{1) factor requires at least a new complex singlet scalar field. If 
the coupling of the new gauge group contains a term proportional to B-L, the absence of 
chiral anomalies demands the presence of additional SM-singlet fermions. The presence of 



right-handed neutrinos, one per generation, removes all anomalies 



Mk2l 



In these minimal models a type-I seesaw mechanism can be introduced, dynamically 



generating neutrino masses. T^ 



with resonant leptogenesis 



DM candidate 



3-50 



le parameters controlling the neutrino masses are compatible 



43 



-|46|. A Z2 symmetry can be introduced to provide a stable 



50|. In these models with type-I seesaw, some fine-tuning is required to 



get two neutrinos almost degenerate and thereby have a successful TeV scale leptogenesis. 
The Z2 symmetry is needed for stabilising the DM, but has no further relation with the 
neutrino masses. 

The inverse seesaw appears more suitable to tackle these issues. Within this framework, 
2 new neutrinos per generations are included. Once the mass matrix (of the left-handed and 
the 2 right-handed neutrinos) is diagonalized, besides the usual light SM-like neutrinos, 3 
pairs of naturally quasi-degenerate heavy neutrinos appear, thereby easily implementing the 



requirements for resonant leptogenesis 51 



51. 



In the context of U {1)b-l extensions of the SM, a model with the inverse seesaw realization 
exists. The model contains two extra dileptons, one of which enters the neutrino mass 

n 

matrix [5^. A Z2 symmetry is advocated that leads to zeroes in the mass matrix. 

A disadvantage of this model is that this Z2 symmetry has to be broken, in order to 
avoid exactly zero-mass neutrinos. The breaking mechanism is not present in the tree-level 
Lagrangian, coming from non-renormalizable operators. 



In this paper we discuss a consistent extension of the SM with a U{1) gauge group, re 



to B-L number, providing for an inverse seesaw mechanism. In contrast to Ref. 5J], we 



ated 



only use renormalizable operators. The mechanism is natural in the sense that we allow for 



all renormalizable terms in the Lagrangian, consistent with the symmetries of the fields. We 
add pairs of fermions with fractional lepton number, so-called "leptinos" to the Lagrangian. 
One of them is odd under an additional Z2 symmetry, the other even as are all ordinary 
SM particles. Both fermions are needed in order to cancel anomalies. The Zi symmetry, 
together with the B-L charge assignments, restricts the form of the neutrino mass matrix. 
At the same time it stabilises the odd leptino, that becomes the dark matter candidate in 
the model. 

We will present in detail a version with only one extra pair of leptinos, which is sufficient 
for the discussion of dark matter. The possibility of successful leptogenesis requires the 
extension of the inverse seesaw mechanism to more generations, in order to provide the 
necessary large phases driving CP-violation. The detailed study is beyond the scope of this 
paper. However, we make some comments in the last section. 

The paper is structured as follows. In the next section, the model is presented. Section UTTl 
collects results for the dark matter abundance generated by the model. In section |IV] the 
possibility of a successful leptogenesis is outlined. Finally, in section |V] we present our 
conclusions. 



II. THE MODEL 



We base our extension of the SM on the minimal Z' model 29|-|33| . The SM gauge 
group is extended by including a f/(l) factor, related to the B-L number, with generic 
mixing with the ?7(l)y. A SM singlet complex scalar x is required for the spontaneous 
symmetry breaking of the further f/(l) group, thereby providing the Z' boson a mass. 
The requirement of anomaly cancellation is fulfilled by introducing one right-handed (RH) 
neutrino per generation. Furthermore, the inverse seesaw mechanism needs extra SM singlet 
fermions, coming in pairs in order not to spoil the anomaly cancellation. Minimally, just 
one extra pair of fermions is sufficient to provide a DM candidate. 

The classical gauge invariant Lagrangian, obeying the 5't/(3)c x SU{2)l x f/(l)y x 
U{X)b~l gauge symmetry, can be decomposed as: 

=Sf = ^YM + =S?s + -Sf/ + -Sfy • (1) 
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A. Gauge sector 



In the gauge field basis in which the kinetic terms in ^ym are diagonal 33|, the covariant 
derivative reads: 



D, = d, + igsT'^G^^ + igwT'^W;' + igiYB^ + i^g^Y + gBLYB-L)B'^ . 



(2) 



This generic model describes a continuous set of minimal UJX) extensions of the SM, that 



can be labelled by the charge assignments of the particles [33|] . As any other parameter in 
the Lagrangian, g2 and gsL are running parameters, therefore their values have to be set 
at some scale. Special sets of popular Z' models (see, e.g. [29|, l55|) can be recovered by 



imposing relations between g2 and gsL- However such relations can be changed through the 
renormalization group running of the coupling constants. 



It is instructive to consider the renorma^ 
couplings. The one-loop RGEs read 



28 



31 



ization group equations (RGEs) for the gauge 



56| 



d 

dt'' 



d_ 
di 



gsL 



d 



1 

167r2 
1 



:92 



[A^^gl, + 2A^^gl,g2 + A^^g^.gj] , 

[A^^g^ {gl + 2gl) + 2A^^ g^dgl + gl) + A^^ gl.g,] 



(3) 
(4) 

For the model we are discussing {Y is the SM weak hypercharge, X = B-L is the B-L 
number), the coefficients are: 



A 



YY 



41 



A 



XX 



32 + (K 



B-L\2 
X ' 



-Nf 



A 



YX 



16 

y 



3 27"^' 3 ' 

if generations of leptinos (section C) are included. Notice the small difference with respect 
to Ref. 311 in the Al^^ coefficient, due to the SM Higgs in the counting. The equations can 
be solved algebraically 3l| : 

1 



2y4^^t = constant , 
2Ay^g, + 2Ay^gsL 



9l 



9i9bl 



A''^g! + gl)-A^^gl, 



constant 



constant 



(7) 
(8) 
(9) 



Particularly interesting is eq. ([H]), which leads to an infrared (IR) fixed point for the model 

Alg2 + ?>2gBL = . (10) 



This fixed point is independent of tlie additional matter we consider in tfie model. The 
reason for this is that the extra fields we introduce (RH neutrinos, leptinos and the singlet 
scalar) are all singlets under the SM gauge group, hence not entering in the diagonal and 
mixed hypercharge coefficients and A^^, that are the only terms appearing in eq. 
In other words, this IR fixed point is a model independent property of the minimal Z' model. 

For the following study, it is more important to focus on measurable observables, one of 
which being the Z' total width. In the approximation of massless fermions and neutrinos, 
with Ng generations of leptinos, of which just the CP-odd leptinos (S'2, the DM candidate) 
are massive, with in first approximation degenerate masses Mgj, the total width reads 




127r \ \ 2 18 V V Mz' J j \ 8 y V 2 2 ^ 

(11) 

Since the S'2 particle is only right-handed, its hypercharge is zero (see table [Tll), so that 
its coupling to the Z' boson does not depend on g2. Moreover, in the evaluation of the 
DM relic abundance with S'2 as the DM candidate, we will see that a resonant annihilation 
with the Z' boson is required. In these conditions, 2M52 ~ Mz' and BR(Z' — )■ S2S2) — )■ 0. 
Being around the resonance, the only parameter that can influence the S2S2 -> Z' partial 
amplitude, and thus the relic density, is therefore the Z' boson width. The smaller the total 
width, the higher the relative partial amplitude. Hence, the minimum relic density, when 
one keeps Mz' fixed, can be obtained by minimizing the Z' width of eq. (ITT]) with respect 
to g2- By direct computation, we obtain that the total Z' decay width is minimized for 

gT''{9BL) = -2^^^gBL. (12) 

When A''^ = 1 (the case discussed here), (7^*" = — 10/13 (^bl- Other important cases are 
for = 3 (discussed in section HVl in connection with leptogenesis) , for which we obtain 
grnin _ —2Q/37 qbl, and A''^ = 0, where the minimum width is for the SO(10)-inspired U{1)^ 
model, (7™" = —4:/5gBL- 

So far, we have neglected that the Z' boson is in general mixed with the SM Z boson. 



Typical bounds from 
be less than O{10-^) 



EP-I measurements at the Z-boson peak require the mixing angle to 



3- 



Further, from a combination of LEP-I and LEP-II data the ratio mass-over-coupling is 



bounded to be bigger than several TeV. In Ref. 



32| the authors reanalyzed the LEP data 



for a model with the same gauge sector as ours, while a specific bound for g2 = can also 



be found in Refs. 



29 



58 



Before moving on to the results, we present here (see figured]) the 95% C.L. exclusions at 
the LHC in the gBL — 92 plane, based on the CMS data at ^/s = 7 TeV for the combination of 
4.7(4.9) fb^^ in the electron (muon) channels 59[. ATLAS has as well published an analysis 
for ~ 5 fb~^ 60|, but their limits are less tight than the CMS ones. Therefore, we will not 



present them here. 
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FIG. 1: Z' exclusions from CMS data, at ^/s = 7 TeV for the combination of 4.7 fb~^ in 
the electron channel and 4.9 fb"^ in the muon channel. The dotted black lines refer to the 

two main benchmark models of this analysis. 

Table H] collects the maximum allowed qbl coupling per given Z' boson mass for the B-L 
model ((72 = 0) and for the ^f™" case of eq. (fT2|) . for A''^ = 1. 

B. Scalar sector 

In order to give the Z' boson a mass through spontanous symmetry breaking, an extra 
complex scalar x is introduced. The scalar Lagrangian reads 

= {D'^Hy D,H + {D^x)^ D,x - V{H, x) , (13) 
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Mz' (TeV) 


gBL for 


gBL for B-L 


2.5 


> 0.8 


0.75 


2.2 


0.81 


0.45 


2.0 


0.58 


0.31 


1.8 


0.39 


0.22 


1.5 


0.24 


0.13 


1.0 


0.13 


0.08 


0.5 


0.05 


0.03 



TABLE I: 95% CL. exclusions for the benchmark models of interest. Couplings smaller 

than those in the table are allowed. 



with the scalar potential given by 

ViH, x) = m'H^H + fi'\x\' +)^i{H^Hf + | x 1' +^:^H^H \x\\ (14) 

where H and x ^^e the complex scalar Higgs doublet and singlet fields, respectively. 

The scalar sector is now made of two real CP-even scalars, a light one hi and a heavy 
one /i2, that are mixtures of the Higgs doublet field and the singlet field: 



^h^ 


H 







cos a — sm a 
sin a cos a 




(15) 



where v{v') is the VEV of the doublet (singlet) field, while the mixing angle — f < a < | 
can be expressed as: 

tan 2a = ^^^^ , (16) 

with Ai 2,3 being the parameters entering in the quartic pieces of the scalar potential. The 
presence of extra heavy neutrinos and the Z' boson will alter the properties of the Higgs 
bosons. The scalar mixing angle a is a free parameter of the model, and the light (heavy) 
Higgs boson couples to the new matter content proportionally to sina(cosa), i.e. with the 
complimentary angle with respect to the interactions with the SM particles. 

The B-L -breaking vev v' can be expressed in terms of the Z' mass and the gauge couplings 
as follows: 
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Here we have chosen the B-L charge of the complex singlet scalar to be Y^~^ = |, as 
determined by our choice of the Yukawa sector. 

The LEP bounds on Mz' / Qbl can provide an absolute lower bound on the VEV v' through 
eq. ( ITTj) . For the benchmark scenarios of interest here, the bounds are: 

^2 = 0: t/ > 9 TeV , (18) 
t;'>6.7TeV. (19) 



C. Yukawa sector 



Moving to the Yukawa sector, the core of this paper, the inverse seesaw mechanism can 
be implemented by means of i extra pairs of RH fields beside the usual 3 RH neutrinos {vr), 
the latter required by the anomaly cancellation. With only renormalizable operators, the 
Yukawa Lagrangian for the neutrino sector reads 



C'y = -y^iLHvR - y^'v'kSlx - yf{s^'S{x^ - yf{s\YSix + h.c. , (20) 

where 5*2 are the only odd fields under a Z2 symmetry that is introduced to avoid unwanted 
S1-S2 mixing. A summation over z, j is implied and an identical number of Si and S2 fields 
is required by anomaly cancellation. As a concrete model for the DM study, we focus on 
i = j = 1 case. By convention, we choose to implement the inverse seesaw mechanism in 
the third generation of leptons only. In the following, the formulas will refer to the latter 
and we drop the "z" index, unless otherwise specified. 

The gauge invariance of the Yukawa Lagrangian is achieved by solving Y^'^ + Y^ '^ = 
-Y^^^ and 2Y^-^ = ±Y^-^, with Y^-^ the B-L charge of the fields. The following two sets 
of solutions exist, considering that Y^~^ = 1 is fixed by anomaly cancellations: 



i) Y<^^~^ = —1 and Y^~^ = 2, that is a replica of the type-I B-L model '< 



. In this case. 



Majorana masses for all the neutrinos would be allowed, as well as a iiHSi term; 

ii) Y^^~^ = 1/3 and Y^~^ = 2/3. This is the new solution proposed here, that forbids all 
unwanted terms. Because the B-L charge is one third of that of the normal leptons 
and because they do not carry a colour charge, the fields will be called "leptinos" 
throughout the rest of this paper. 
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^2 
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SU{2)l 
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B-L 


-1 


-1 


-1 
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1 
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2 
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Z2 
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+ 
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+ 
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+ 



TABLE II: Y and B-L quantum number and Z2 parity assignations to chiral fermion and 

scalar fields. 

Table |IT1 summarises the particle content and the charge assignments. 

After spontaneous symmetry breaking, the two Higgs fields acquire vacuum expectation 
values (called v and v', respectively). As a consequence, the first 3 terms of eq (120!) lead to 
the following mass matrix for the interacting neutrinos, in the {uf, uji, Si) basis: 



^ Md ^ 

Ml Mn 
y Mn Ms, J 



(21) 



where 

Md = ^v, Mn = ^v', Ms, = V2y'\/. (22) 

When all 3 generations of leptinos are considered, y'^ and y^^ are in general 3x3 complex 
matrices, and it is not restrictive to consider as a 3 x 3 real and diagonal matrix. 
After diagonalization of eq. f l2T]) . the 3x3 neutrino mass matrices are 

M,^ ~ MnM^'Ms,{M^)-\Mny , (23) 
Ml = Ml ^Ml + Ml, (24) 

However, it is beyond the scope of this paper to consider the most general case in detail. 
We limit ourselves to the model with only one generation of leptinos, which is largely suffi- 
cient for the dark matter question. In this case, one can without loss of generality choose a 
basis in which both the 3-component vector y^ and the single parameter y^^ are real. The 
above formulas simplify, since the seesaw mechanism acts only on the third generation, so 
that is a 3 X 3 matrix and eq. (123|) gives the masses of the 3 neutrino eigenstates, where 

and combine in a quasi-Dirac fermion. Regarding the first and second generations, 

10 



the LH and RH neutrinos obtain the usual Dirac mass term, for which (9(10~^^) Yukawa 
parameters are required to fit the hght neutrino masses. 

The mass matrix of eq. (^T^ is diagonahzed by a 3 x 3 unitary matrix. Such a matrix can 
be parametrized by 3 angles and some phases. The latter will be neglected here. Formally, 
the following generic parametrization can be employed: 





fl 
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Cl3 





'S'la 


\ 




C12 


'S'12 


o\ 




u = 







5*23 









1 









— 5*12 


C12 





(25) 




1° 




C23 
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— Si3 





Cl3 
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V 





where Sij{Cij) = sinajj(cosajj). In good approximation, the angle defining the mixing 
between and Si is found to be very close to maximal, i.e., a23 ~ vr/4, while the angle 
parametrizing the mixing between z/^ and 5*1 is very small: S'la ~ MijMs^/Mff ^ 1. The 
unitary matrix of eq. fl25l) can then be simplified to 







C12 


S12 


o\ 


u - 




5l2 

V2 


1 C12 


1 




V 


S12 
V2 


C12 

V2 


72 / 



(26) 



where S'12 ~ = controls the mixing between ul and uji. Altogether: 

Mat yiyv' 









\ 




= U 










^ S'l 


/ 



(27) 



Neglecting intergenerational mixing, eq. (123|) can be rewritten as 

2 2 

M,, ~ MDM^''MsAMl)-\M*^f = V2^^^ . (28) 

Vn ^ 

To obtain light neutrino masses compatible with experiments (i.e., My^ < 1 eV), considering 
v' = 10 TeV, eq. simplifies to 

y.if^V<10-^°. (29) 

A possible solution, which allows for the heavy neutrinos in the O(IOO) GeV range (suitable 
for the LHC), requires vn ~ O{10-^). Then, ~ 0(10-^) implies y^i < O{10-^). 
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A fundamental difference with the model of Ref. [5J] exists in the non-interacting neutrino 
sector, the 5*2 fields (here we consider only one of them). In eq. f l20|) it is clear that the 5*2 
field acquires a mass after the U{1)b-l symmetry breaking when x gets a vev (called v'): 

Ms, = V2y,2v'. (30) 

Therefore, its mass is expected to be of (9(1) TeV, if also v' is at TeV scale. The Z2 symmetry 
protects this particle from decaying, making it a suitable dark matter candidate. 

Given that the 5*2 field acquires a mass through a Yukawa term similar to the Majorana 
mass of the Si field (see eq. ( 130|) ). it is expected that their Yukawas {ysi and ys2) lie in 
similar ranges. As we will see, the 5*2 field is a suitable DM candidate if its mass is of the 
order of hundreds of GeV, that requires ys2 ~ O(10~^), similar to the order of magnitude 
for required from light neutrino mass constraints, as one would naively expect. 

Such range in masses allows for a suitable relic density due to resonant decays via the 
heavy Higgs boson only. Relaxing the latter requirement ends up in a larger variety of phe- 
nomenological implications, such as heavier DM candidates and resonant annihilations with 
the Z' boson, as well as more interesting LHC phenomenology in the interacting neutrinos 
due to a larger mixing with the LH fields, because of a larger coupling of (9(10~^). 



III. RESULTS: DARK MATTER 

The 5*2 field as defined in the previous section is a suitable candidate for WlMP-like 
cold dark matter because it is odd with respect to the Z2 symmetry and thus stable on 
cosmological time scales. An abundance of S2 particles which is thermally produced in the 
early universe will survive after freeze-out from the thermal bath until today to make up 
the observed dark matter in the universe. 

To compute the relic density of S2 dark matter we used the program MicrOMEGAs 



62| . in which the model has been implemented via LanHEP |63|. The remaining numerical 



61 



analysis has been performed in CalcHEP |6J]. The a priori unknown Yukawa couplings of 



the extra neutrino fields turn out to be negligible in the calculation of the relic density. For 
concreteness, in the following we fixed = 10~^, = 0.06 and y'^^ = 10~^, owing to 
m,^i = 1.5 ■ lO^^'' and M^^ = M^i^ = 636 GeV for f' = 15 TeV. Further, the Higgs mixing 
angle has been set to sin a = 0.1. 
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FIG. 2: DM relic abundance as a function of the candidate mass for M^^ = 800 GeV and 
Mz' = 2 TeV, for two different choices of the scalar mixing angle. For the other 

parameters, see the text. 



The recent discovery at the LHC 65| if interpreted as a SM-like Higgs boson, does not 
constrain this parameter very strongly yet, but implies it is small. The neutrino couplings 
and scalar mixing angle a will stay fixed to these values throughout the following analysis. 
Therefor our free parameters are Mz', Mh,^, Qbli 92 and the dark matter mass Ms2- 

For a first assessment of the situation we fix qbl = 0.1, g2 = 0, = 0.8 TeV and 

Mz' = 2.0 TeV. We expect the coupling qbl and the Higgs mixing angle a to have an 
impact on the relic density, as they control the interaction strength of the Z' boson and the 
Higgs bosons, respectively, with the DM particle and the remaining particles of the model. 

We compute the relic density as a function of the 5*2 mass, as shown in figure |2j Comparing 
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the resulting graph with the observed abundance of dark matter: 



QnMh' = OAn7t'oS, (31) 

it turns out that for ^ Mz'/2 and Ms^ ^ M//2/2, the rehc density is always too large 
by several orders of magnitude. There are two ways to reduce Qh"^: one is to lower the 
vev v' (which results in an increase in the Yukawa couplings), which will make the heavy 
Higgs resonant annihilation more efficient. The other possibility is to increase the coupling 
qbl to enhance the Z' annihilation efficiency. However, an arbitrarily small relic density is 
forbidden, due to the LEP lower limit on v' of eqs. (fT8l) - (fT9l) . and by the upper limit on 



66|. Obviously, the 5*2 Yukawa coupling is hnked to the v' and 



qbl due to RGE analysis 
therefore to qbl- 

We conclude from this first assessment that for the 5*2 particle to be the dark matter, it 
must annihilate efficiently via a resonant heavy Higgs or via the Z' boson. Also the precise 
value of a is not relevant in the resonant regions, that are the only important regions for 
our analysis, as long as it is small, as expected if one relies on the confirmation of the LHC 
hints on a light SM-like Higgs boson. In the following we will study these two different 
mechanisms, fixing sin a = 0.1. 



A. Higgs boson resonant annihilation 

For small DM masses, where the main annihilation channel is via the heavy Higgs boson, 
the relic density is proportional to the 5-L-breaking vev, since Qh'^ oc - oc oc v'"^, via 
eq. (|30l) . According to LEP, the vev is constrained from below as in eqs. ( fT8l) -( !T9l) for the 
models of interest. The left panel of figure [3] shows that for Ms^ — ^^2/2 the minimal 
abundance just matches the observed value. 

In the right panel of figure |3] we see that the demand, that the resonant S2 annihilation 
proceeds via the Higgs channel, strongly constrains from above the vev v' as a function of 
the heavy Higgs mass M/^-j, v' < 11.7 TeV. It is interesting to note that this upper bound 
on v' is independent of the scalar mixing angle a. At the same time, the heavy Higgs mass 
is constrained from above depending on the value of a. For instance for sin a = 0.1, < 
1.8(3.0) TeV for ^2 = 0(^r")- 
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FIG. 3: (Left) DM relic abundance as a function of the DM mass for = 800 GeV, for 
some values of the vev v' and sin a = 0.1. The blue shading represents values of the vev for 
which an allowed DM mass exists, for this Mh2, as taken from the right panel. 
{Right) Allowed region for v' for which a DM mass yields the correct relic density, as a 
function of the heavy Higgs boson mass, for the three values of sin a = 0.05, 0.1, 0.3. The 

LEP exclusion is as in eqs. (fT8l)-( fT9|l . 

B. Z' boson resonant annihilation 

Due to LHC direct searches, the Z' boson mass has to be above 2^2.5 TeV, depending on 
the gauge mixing coupling g2. For the Z' annihilation mechanism to be effective, as already 
stated, a resonant condition has to be matched, meaning that the DM candidate has to be 
rather heavy. However, the resonance decay can still be very effective. 

Figure H] shows in more detail that only around the Z' resonance the annihilation is 
sufficient to match the abundance constraint. For simplicity, we fix here g2 = 0, so that the 
relic density is inversely proportional to the square of the gauge coupling, i.e., Qh"^ oc (~^- 
The lower limit on the vev v' translates into an upper limit on the coupling, for a fixed Z' 
mass: in figure Hlthe Z' mass is 2 TeV, which gives the upper limit gsL < 1/3, via eq. (ITTI) . 
Then, the demand Qh"^ ~ 0.1 gives a lower limit for the gauge coupling, gsL > 0.21. 

Although the gauge coupling gsL is expected to give the major contribution to the relic 
abundance evaluation, the impact of g2 might be not negligible either. For fixed Mz' and gsL 
values, we study here the effect of the mixing gauge coupling. We choose two different sets 
of Mz' and gsi such that for (72 = the relic density just satisfies the abundance constraint. 
Figure [5] then shows that g2 can have an impact on the relic abundance, especially it can 

15 




10000 F 




0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 



Ms, [TeV] 

FIG. 4: Relic density as a function of the DM candidate mass around the Z' peak. The 
lower curve, for qbl = 0.5, is for illustrative purposes only, being already excluded by LEP. 
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FIG. 5: Variation of the relic density (at the Z' resonance) with (72 for a choice of qbl (see 
the text for further details), for (left) Mz' = 2 TeV and (nght) Mz' = 3 TeV. 

lower the latter when assuming small negative values. 

A general feature is that 0,{g2) is growing for \g2\ — )■ 1, so that there exists a value of 
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FIG. 6: Existence of a suitable DM candidate: the allowed region is the one above the 
dashed curves, in the qbl — Mz' plane, (left) for g2 = and (right) for g2 = gT^^'igBi), 

which minimizes the Z' width, hence allowing for a smaller gsL (and therefore a smaller 
cross section) per fixed Mz'- The red shading combinations are forbidden by LEP 

(eqs. ( !T8|) -( !T9|) . respectively), the black (solid) lines are the LHC exclusion, as in table E 



g2 that minimizes ^l{g2)- This value of g2 corresponds to the (7™" of eq. ( 1T2|) . that also 
minimizes the Z' width. Fixing g2 to g^^^{gBL) returns a minimum value for gBL such that 
the abundance constraint can be matched, which is roughly 5% lower than the lower limit 
on gBL in the case g2 = 0. We can now study the relic density as a function of the Z' boson 
mass. This is done in figure El that shows the range of allowed values for gBi as a function of 
Mz'- The curves, for both g2 = and g2 = g2^^"{gBL), hmit the existence of a DM candidate 
mass with suitable relic density above the curves themselves. 



As a result, the Z' mass is constrained to Mz' < 3.5(5.0) TeV for the (72 = and the 
g2 = gT^"'{gBL) case, respectively, which also imposes an upper bound on the gauge couplings 
when combined with the exclusion limits from LEP. Notice that this value for the Z' mass 
is within the LHC ultimate reach 4l|. Since we must be on a resonant regime, the upper 



bound on the Z' mass translates on a upper bound for the 5*2 mass, M^^ < 2.5 TeV. 
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FIG. 7: Spin-independent direct searches for maximum Yukawa couplings ~ ^s^l'^' ■ 
Only allowed masses are plotted, from figure right), for = g^''"'- 



C. Direct detection 



The spin independent leptino-nucleon interaction is mediated by the light Higgs boson 
and is therefore sensitive to the value of (sin a cos a)^. Figure [7| shows that even for the rather 
large value of sin a = 0.3, the leptino-nucleon cross section is several orders of magnitude 



below the actual exclusion limits from XENONIOO [67|, at most of roughly 10~^^ cm^ for 
a DM candidate mass of around 50 GeV. Since we know that sin a has to be small, also 
future direct detection experiments will not be able to restrict or to detect the leptino. We 
checked that the exchange of Z' bosons consistent with LEP exclusion limits gives rise to 
smaller cross sections than Higgs exchange, hence we did not further consider this process. 
Even easier to evade are the constraints from spin-dependent leptino-nucleon interaction 
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experiments. Here the only mediator that can play a role is the Z' boson, too heavy for 
these cross sections to be of any interest. 

D. Extension to families 

As we will discuss in the next section, a successful leptogenesis needs > 2 to provide 
the required large CP-violating phases. Here we want to comment on the impact of extra 
leptino generations in the DM analysis carried out so far. The extra neutrinos affect the 
Z' width only, see eq. (fTTj) . and will be discussed later. First we focus on the leptinos. If 
the 5*2 fields possess a mass hierarchy we have a situation similr to SUSY dark matter. The 
heavier particles decouple earlier and decay rapidly to the still thermalized lighter particles. 
An observable effect appears only if the relative mass difference is less than 10%. If they 
are mass degenerate, however, we get fi^v^ = A^^ • fii. Thus to have resonant annihilation 
the couplings have to be increased by a/A^. Figure [8] shows how this affects the parameter 
space. For resonant heavy Higgs annihilation the vev v' must be reduced by a factor a/A^, 
so that only the cases A^^ = 1,2 can still match the relic abundance for Higgs resonant 
annihilation, while A^ = 3 is just touching the LEP exclusion limits (for = (7™") from 
below and therefore such resonant mechanism is not sufficient. 

The case of resonant Z' annihilation works well also for A"^ = 3. However, the allowed 
parameter space becomes tighter for higher A^. These results are valid when all the leptinos 
are mass-degenerate, and represents the worst possible case. All others, i.e. when just 2 
are degenerate or with a tight mass hierarchy, will be somewhere between the case of 1 
generation and the case of 3 generations, degenerate in mass. 

We turn now to study the effect of having 3 generations of leptinos, as required for 
leptogenesis, on the dark matter relic abundance due to the extra heavy neutrinos. The 
model with 3 generations of leptinos has 3 light and 6 heavy Majorana neutrinos in total. 
The one with only 1 generation of leptinos instead accounts for 3 light neutrinos, 2 of 
them being Dirac particles, and 2 heavy Majorana neutrinos. Although the total number 
of relativistic degrees of freedom, commonly addressed as g*, is basically unchanged, the 
proliferation of neutrinos affects the Z' width, and this could impinge on the evaluation of 
the relic abundance at the Z' resonance. For the same setup as in section UTTt figure [9] shows 
the total Z' width in the 2 different cases. 
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to Ni = 3) and (b) percentage variation between the 2 models in the qbl — 92 plane, for 
Mz' = 2.0 TeV. For the heavy neutrino masses, see the text. 



Despite the larger numbers of possible final states into which the Z' boson can decay for 
Ni = 3 with respect to Ni = 1, the latter case has a larger Z' boson partial width into light 
neutrinos since 2 of them are Dirac particles. Moreover, finite mass and threshold effects 
for the heavy neutrinos diminish their branching ratios. Altogether, this induces an almost 
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exact compensation between the larger number of heavy neutrinos for A'^ = 3 and the larger 
partial width due to the 2 Dirac light neutrinos for A^^ = 1. In the chosen setup for the 
Yukawa couplings, the relative variation of the total Z' width is below 1% in most of the 
parameter space, getting above 5% only in a limited region around small values of the gauge 
couplings. In these plots, heavy neutrino masses are rrihy = 637 GeV for qbl = 0.2, and 
scale inversely with the latter. When their mass is above Mz'/I (or slightly before due to 
threshold effects, i.e., here for qbl < 0.15), the Z' — )■ A^A^ (A^ = Uh, v'l^ channels are all 
suppressed or simply forbidden and the above mentioned compensation is not taking place. 
Even though the model with 1 generation of leptinos has here a bigger Z' width due to the 
larger partial widths into the light Dirac neutrinos, the excess is never above 20% of the 
total Z' width in the case of 3 generations, and only in a tiny corner of the parameter space, 
where both gauge couplings are small. 

As intimated, the variation of the total Z' width is very small, mostly below 1%, so that 
the impact on the DM relic abundance is also negligible. This is because, as we can see by 
comparing figure (Hl^a) to figure [S](left) (by reading the value of the total width in the first 
figure and the respective value for the relic abundance in the second figure, for qbl = 0.2 
and Mz' = 2.0 TeV, as a function of g2), the relic abundance scales with the square root 
of the Z' width. A 1% variation in the total width determines a 0.5% variation of the relic 
density. 

We can conclude that the DM study for Ni = 1 of section IIIII is independent of the 
particular value of Ni if the leptinos are not mass degenerate. 



IV. FURTHER RESULTS 

We comment here on further possible implications of the model. First, the possibility of 
leptogenesis is described, necessarily requiring the presence of more than one generation of 
leptinos. Next we comment on how our model has an effect on the effective number of light 
degrees of freedom A^e//, measured in cosmology. Extra degrees of freedom consist of the RH 
neutrinos, forming Dirac fermions with the LH counterparts when less than 3 generations 
of leptinos are present. 
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A. Leptogenesis 



As discussed in the introduction, the inverse seesaw mechanism is a suitable mechanism 
for a large CP asymmetry even at the TeV scale due to the naturally small mass splitting 
between the heavier neutrino eigenstates, a result that is built in into the model and does not 
require fine tuning. Hence, it is a favourable model for the so-called resonant leptogenesis. 
Despite a large loop enhancement due to resonance, large phases are still required in the 
off-diagonal terms oi Ud to obtain an 0{1) CP asymmetry. In the A^^ = 1 model discussed 
so far, the first 2 generations of light neutrinos require (9(10~^^) Yukawas, which in turn 
also means that the CP-violating decays of the heavy Majorana neutrino pair are similarly 
suppressed. Large phases but small masses can be achieved extending the seesaw mechanism 
to the other generations. 

We will show that leptogenesis, compatible with neutrino masses and mixing, is possible 
in the model. The CP asymmetry is generated by the decays of the heavy neutrinos: 

1 Im {yDVDfij 

requiring large phases in the off-diagonal elements of the Dirac Yukawa matrix yo- In our 
model, these large phases, compatible with neutrino data, are possible for Ni > 2. The 
lepton number violating loop factor f^J, when quasi- degenerate heavy neutrino pairs are 
considered, is 

- Mf 

f = i 1 ('i'i) 

{Mf-Miy + {M,T,-M,T,y ^ ^ 

In the inverse seesaw case, Mj ~ Mj and Tj = Fj = F are naturally recovered, so that f^j 
can easily be 0{1). Notice that no fine tuning is required. 

Once an asymmetry in the lepton sector is produced, electroweak sphaleron processes 
take place and move the asymmetry to the baryon sector 
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Vb = J^Vb^l ■ (34) 

Altogether, the final baryon asymmetry can be written as 

?7b ~ 10"^ (2; oo) , (35) 

where the 10~^ pre-factor accounts for the sphaleron efficiency and for photon dilution after 
recombination. The CP asymmetry Si was defined in eq. ( 132|) . while Ki is the efficiency 
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factor obtained after solving the relevant Boltzmann equations: 

Ki{z) ~ / dz' ^'[^ ^ ^}'^'ln\Tenf a X ^Xp 



zO 



dz' D{K„z')+4Sz'K'{z') 



dz"Win{K,,z")5l 



(36) 



with 5i = < 1 suppresses considerably the inverse decay (ID) wash-out. 



r 

All the quantities are defined in Ref. 



53| . Particularly important is the Sz' term, the Z' 



scattering processes, which induces a wash-out of the final asymmetry. It has been verified 
that these Z'-induced wash-out processes do not have a severe impact on the final baryon 
asymmetry. This is due to the nature of the inverse seesaw mechanism, that allows for 
large Yukawa couplings, overcoming the Z' processes. A similar conclusion was reached 



in Ref. 



46| . in which the authors explicitly showed that Z' processes do not spoil the 



leptogenesis if the Dirac Yukawa couplings are sufficiently large, as in the inverse seesaw 
case under examination. 

In conclusion, following the similar case studied in Ref. 53|, the Dirac Yukawa matrix 
yo contains all information required to study the leptogenesis, entering both in Si and in F. 
It is sufficient that some off-diagonal elements in y£, are large and complex for leptogenesis 
to be possible. For definiteness, we have verified that the choice of the matrix in Ref. 53| . 
that for A'^! = 3 is a possibility in our setup, does yield the correct baryon asymmetry in our 
model. We stress again that a similar choice for yo compatible with a successful leptogenesis 
is possible also for Ni = 2, even though only one off-diagonal element will be large in this 
case. This is in fact identical to Ni = 3 in the 1-flavour approximation. 

Hence, we have proven that a successful leptogenesis, compatible with the observed pat- 
tern of neutrino masses and mixing angles is possible in our model when at least two gen- 
erations of leptinos are present. However, the complete analysis of the leptogenesis in our 
model, as for instance the detailed comparison of A^; = 2 and Ni = 3 cases, or the impact of 
flavour effects outside the simple 1-flavour analysis, is beyond the scope of this paper and is 
left for future work. 



B. Impact on N^jj- 

As a last application of our model, we describe the possible implications on A^gjj. This 
observable counts the relativistic energy content in the universe at the time of the last 
scattering surface in terms of an effective number of neutrino species. Indications from 
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cosmology result in the observed value of N^^^^ = 4.56 ± 0.75 when combining WMAP, the 
Atacama Cosmology Telescope, baryonic acoustic oscillations data and the measurement of 



the Hubble parameter Hq [68|. Notice that the SM LH neutrinos only would yield 



eff 

3. This mismatch is typically interpreted as an indication of the existence of some extra 
relativistic degrees of freedom that effectively contribute as one unit of N^jrj. In our setup, 
the only extra degrees of freedom that can be relativistic at the last scattering surface are 
the RH neutrinos, when they have only Dirac mass terms, i.e., for less than 3 generations 
of leptinos, as described in section IIIIDI On the other hand, as we have seen previously, a 
successful leptogenesis requires more than one generation of leptinos. If one would like to 
explain both leptogenesis and AN^jrj in our model, less than three generations of leptinos 
should be considered, given that = 3 leads to an effective number iV^jy ~ 3 as in the 
standard model. 

An estimate of the impact on iVe// done as follows. First, the decoupling temperature 
and the contribution to N^^j^ depend on the parameters mz>, Qbl, 92, given that only the 
Z' boson can keep the RH neutrinos in thermal equilibrium with electrons, as the RH 
neutrinos have vanishing hypercharge. Naively, the bigger the Z' cross section, the lower 
the decoupling temperature, which in turns also means that less degrees of freedom are 
relativistic at decoupling. Overall, this increases the RH neutrino contribution to N^jrj. 
Close to LEP exclusion limits, we get a minimum decoupling temperature of 380 MeV. 
At this temperature the number of relativistic degrees of freedom is 20 if the QCD phase 



transition takes place at 450 MeV 



69l | . This yields the highest value for AN^jjr in our model: 



0.18(3 — A'^). For lower QCD phase transition temperatures, the numbers of relativistic 
degrees of freedom rapidly increases, suppressing the impact on N^jj: to the percent level or 
below. 

Thus we see that although the model can contribute to N^jf, its impact is only marginal, 
especially if one would want to implement a successful leptogenesis. In fact, A^ > 2 is 
required by latter, while A^ — )■ maximises AN^jrj. 



V. CONCLUSIONS 



We constructed a simple model in the class of the minimal Z' models, where the extra 
f/(l) gauge field is coupled only to hypercharge and B-L . We added right-handed neutri- 
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nos for each generation of fermions, as required by the absence of chiral anomahes in the 
theory. Beyond this we added extra pairs of leptons with fractional lepton number, which 
we therefore called leptinos. One of the leptino in the pairs is chosen to be even and the 
other one to be odd under an additional Zi charge. We were able to construct an inverse 
seesaw mechanism for neutrino masses. The mechanism is natural in the sense that all 
possible terms consistent with the symmetries of the theory are present in the Lagrangian. 
Nonetheless some of the entries in the neutrino mass matrix are zero, because of the choice 
of the representations. The reason is the presence of a fractional lepton charge, which is the 
new feature of the model. 

The odd leptino is a candidate for dark matter, as it is weakly interacting, massive and 
stable. We have shown that the correct dark matter density can be generated if the leptino 
is annihilated through a resonance by either the Z'-boson or the Higgs-boson, related to the 
breaking of the B-L symmetry. We studied limits on the parameters of the theory, coming 
from LEP, hadron colliders and the dark matter abundance. We found that the limits are 
such, that the Z'-boson and extra Higgs-boson lie within the range of the LHC, though the 
full design luminosity might be needed. The cross sections are too small for present direct 
search experiments for dark matter. Resonant leptogenesis is possible in the presence of 
more than one pair of leptinos, but contains no particularly new features compared to other 
models of leptogenesis. 

In conclusion the model provides a very simple extension of the standard model, contain- 
ing a number of desirable features like dark matter, leptogenesis and (inverse) seesaw. At 
the same time, being a singlet extension, it does not lead to phenomenological problems, 
such as flavour changing neutral currents. 

The model we constructed appears to be able to give a realistic description of cosmological 
data and can be tested at the LHC. However, it is not possible to describe at the same time 
leptogenesis and ^^f^- 
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